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Abstract. The dynamics of motor protein molecules that have two subunits is 
investigated using simple discrete stochastic models. Exact steady-state analytical 
expressions are obtained for velocities and dispersions for any number of intermediate 
states and conformations between the corresponding binding states of proteins. These 
models enabled a detailed description and comparison of two possible mechanisms of 
the motion of motor proteins along the linear tracks: the hand-over-hand mechanism 
when the motion of subunits alternate, and the inchworm mechanism when one 
subunit is always trailing another one. It is shown that particles in the hand-over- 
hand mechanism move faster and fluctuate more than the molecules in the inchworm 
mechanism. The effect of external forces on dynamic properties of motor proteins 
is discussed. Finally, a method is proposed for distinguishing between these two 
mechanisms based on experimental observations. 
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1. Introduction 

Motor proteins, also called molecular motors, are processive enzyme molecules that play 
a fundamental role in most biological processes, but especially in cellular transport, 
motility, cell division, and transcription pfl El Ej. Motor proteins, such as kinesins, 
dyneins, myosins, polymerases, helicases, etc., function by stepping between equally 
spaced binding sites along the rigid polar linear tracks (microtubules, actin filaments, 
DNA molecules) , and the motion is powered by the hydrolysis of adenosine triphosphate 
(ATP) or related compounds. The mechanisms of the transformation of chemical energy 
of hydrolysis into the mechanical work in motor proteins are not yet fully understood 

m 

Recent experimental advances have allowed for the determination of structural and 
dynamic properties of motor proteins with a high-degree precision at a single-molecule 

level H0ini00l3[ini[m[I2[I3d[H[in]- Cr y stal structures suggest that 
many motor proteins have two domains elastically coupled together, each capable of 
hydrolyzing ATP and moving along the linear track [3J E] . Two possible mechanisms of 
the coordinated motion of protein molecules with two motor heads have been proposed 

In a hand-over-hand mechanism, at each step only one motor head undergoes a 
sequence of mechanochemical transitions so that the motor subunits alternate between 
trailing and leading positions at the beginning of the cycle. In this mechanism the motor 
subunits are fully equivalent to each other. In contrast, according to an inchworm 
mechanism, one motor domain is always ahead of another one during the cycle, i.e., 
the motor heads are not equivalent at all times. Experiments on single-molecule 
fluorescently labeled myosins V, that step along the actin filaments, and on kinesins, 
that move along microtubules, support the hand-over-hand mechanism for these motor 
proteins ^21 EH E3 EE] • However, there are indications that dyneins probably utilize 
the inchworm mechanism [T7j . 

Successes in experimental studies strongly stimulated many theoretical investiga- 
tions of mechanisms and dynamics of molecular motors [TBI EH EH EJJ EH Ell EH E3 
I2H1 EH EH EH ED] • Most theoretical work on motor proteins follows two main directions. 
One approach utilizes the concept of thermal ratchets [TH EH ED] • According to this 
idea, the motor protein is viewed as a Brownian particle that moves in two different pe- 
riodic but asymmetric potentials, switching stochastically between them. This method 
takes into account the internal structure and interactions between different domains in 
protein molecules, however, the results are mainly numerical and depend on the specific 
potentials used in calculations. An alternative approach is based on multi-state discrete 
stochastic (chemical kinetic) models HHEDlEUEHEHlElESlEHlEn!. In this method 
the molecular motors are associated with particles that move along one- dimensional 
periodic lattices with different forward and backward rates. The lattices corresponds to 
biochemical pathways for the motor proteins, while the sites in the period describe the 
biochemical cycle when the protein molecule travels between two consecutive binding 
sites. Using this mapping of the motion of a random walker and applying the method 
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of Derrida > exact analytical expressions for the mean velocities and dispersions are 
derived for any number of intermediate states (i.e., for the period of any size) and for 
different complexity of biochemical pathways |2B 122 EH EH ESj • It was demonstrated 
that this approach allowed for the successful analysis of the full dynamics of single ki- 
nesin and myosin V molecules |27)| EH]- However, the weakness of this method is the 
fact that the internal structure of motor proteins, namely, the motion and interactions 
of different subunits, is not taken into account. 

Determining how the different motor heads move relative to each other is critical 
for the overall understanding of motor protein's dynamics and functions. Current 
experimental methods with single-molecule fluorescent labels, that distinguish between 
the different types of molecular motion, require a detailed knowledge of the protein 
structure which is not always available. In addition, the labeled proteins may change 
their biochemical properties in comparison with the original species. However, it would 
be more advantageous to use simpler less-invasive experimental methods along with 
better theoretical models to study the specific mechanisms of molecular motors. In 
this article we investigate the dynamics of motor proteins by developing a set of 
simple multi-state discrete stochastic models. In our approach the motor proteins 
consist of two interacting particles that correspond to different motor subunits in real 
enzymes. Explicit formulas for the velocities and dispersions are obtained for two 
different mechanisms of motion. We suggest a method to distinguish between two 
possible mechanisms by analyzing time trajectories of single motor proteins obtained in 
optical-trap experiments [HI IH1 El EH UJ] in combination with the bulk biochemical 
kinetic data. 

2. Theoretical Approach 

2.1. General Model 

Consider a motor protein molecule with two subunits that travels along the filament 
track. We model this system as two identical interacting particles moving on a periodic 
one-dimensional lattice, as shown in Fig. 1. There are N intermediate discrete states 
on a biochemical pathway between two consecutive binding sites. In the simplest 
approximation, we assume that the particles interact through hard-core exclusion, i.e., 
they cannot occupy the same site. Also, the particles cannot run away from each other. 
If X\ and X2 are the positions of the motor particles then 

\%i — Xi\ < m, (1) 

where m is an integer that specifies how far apart two motor domains can be found 
in the protein molecule. In the lattice the sites x = ±Nl (7 = 0, 1, • ■ •) correspond to 
binding sites of the molecular motor. The distance between two consecutive binding 
sites is d, which is equal to 8.2 nm for kinesins and dyneins moving on microtubules and 
36 nm for myosins V and VI traveling along the actin filaments P3EHH!- 
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Figure 1. General schematic view of periodic stochastic models of motor proteins 
consisting of two subunits. Two parts of the molecular motor cannot occupy the same 
site and cannot be more than m sites apart. The motor domain at site j can make a 
forward or backward steps with the rate Uj or Wj , correspondingly, if these transitions 
are allowed by another motor domain. 



The particle at site j moves forward (backward) with the rate Uj (wj) if the site 
j + 1 (j — 1) is available and the move does not violate the condition (JI}: see Fig. 1. 
Because of periodicity the transition rates are related, Uj±Ni = Uj and Wj±m = Wj for 
/ = 0, 1 • ■ ■ and < j < N — 1. The dynamic properties of motor proteins are specified 
by the drift velocity 

V = V({u j ,w j })= lim (2) 

t^oo OA, 

and dispersion (or diffusion constant) 

D = D({ Uj , Wj }) = hm 1| [<* 2 (t)> - (x(t)) 2 ] , (3) 

where x(t) is the position of the center of mass of the protein molecule at time t. It is 
convenient to express the degree of fluctuations of the molecular motor in terms of a 
dimensionless function called randomness [1] 

r = w (4) 

This function sets bounds on the number of rate-limiting biochemical transitions 
and thus yields an important information about the mechanism of a motor protein's 
processivity 123 H3 • 

The motor proteins in experiments and in the cellular environment frequently work 
against external loads PJ El E]- External forces modify the transition rates in the 
following way El 122 121 IM 123 

u 3 (F) = u J (0)expl-^- F j, ^(F)=^(0)expl+-^j, (5) 

where ^2^=^(0^ + 9j) = 1, and dp are load-distribution factors that specify how the 
external load changes the energy activation barriers for the biochemical transitions from 
the state j. 
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Figure 2. The general picture of the hand-over-hand mechanism for motor proteins 
with two heads, (a) The black particle is a leading motor domain and it cannot move 
forward. The grey subunit is passing over the black subunit. (b) The black particle 
is now the trailing subunit and it is available for the forward motion. Arrows indicate 
the allowed transitions. 



The dynamic properties of the motor proteins with two domains depend on the 
specific mechanism of the motion. Below we consider in detail the hand-over-hand and 
the inchworm mechanisms. 

2.2. Hand-Over-Hand Mechanism 

In this mechanism, the trailing subunit makes N intermediate steps and becomes the 
leading particle, as shown in Fig. 2. Then the next particle makes N transitions. 
During the cycle each head advances the distance 2d so the center of mass of the protein 
molecule moves only the distance d. This mechanism then can be viewed as a motion of 
two particles on periodic parallel one-dimensional lattices, where the distance between 
neighboring binding sites is 2d. Because the particles are identical, this picture is easily 
mapped into the motion of the single particle (center of mass) on the original one- 
dimensional lattice for which the dynamics is well understood p?T| I21j. 

At large times the exact expression for the drift velocity is given by [HH EI] 




(6) 



and 




(7) 



The corresponding expression for dispersion can be written as [KB 1^ 

d fdU N + VS N (N + 2)V\ 




(8) 



where the auxiliary functions are 




(9) 



j=0 



k=0 



3=0 
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Figure 3. Two possible configurations for the inchworm mechanism of the motion with 
m = 2. The black particle is always leading. The allowed transitions are indicated by 



arrows. 



and 



N-lj-k 



■ 1:11 



W i+ i 



(10) 



k=X j-1 



It can also be demonstrated that in this mechanism Th h > 1/N for any set of transitions 
rates {u^Wj} 123 HI]- 

2.3. Inchworm Mechanism 

Now consider the inchworm mechanism of the motion of the motor proteins with two 
subunits. In this mechanism one particle is always leading the other one, as shown 
in Fig. 3. In the simplest approximation, we discuss only the case m = 2 in Eq. 

i.e., the two motor particles can only be found on two nearest-neighbor sites or 
next-nearest-neighbor sites: see Fig. 3. 

To determine the dynamic properties of motor proteins in this model we develop a 
method that generalizes the original Derrida's approach (231221 El]- The first step is 
to introduce Pj t k{l,t) which is the probability to find the trailing subunit of the motor 
protein molecule at state j and the leading subunit at state k (k = j + 1 or j + 2) at site 
/ at time t (see Fig. 3). The time evolution of this probability is governed by Master 
equations: 



dPjj+i(l, t) 
Jt 



Uj-xPj-\j + i{l,t) + W j+2 Pj,j + 2(l,t) 

- (u j+1 + Wj) P jtj+1 {l,t) 

U j +1^3,3 +l^i t) + w j+lPj+l,j+2{h t) 



dt 





Because of the conservation of probability, we require 



+oo N—l 




for all t. 



(12) 
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Next following Derrida's method we define auxiliary functions 



+00 



+00 



(13) 



and 



i=— 00 



+00 



i=— 00 



l=—oo 
+00 



(14) 



i=— 00 



Using the Master equations (JTTJ) we derive 

^Bjj+iit) = + w j+2 B jJ+2 (t) - (u j+ i + Wj) B jd+1 (t), 

d 



-^B jd+2 (t) = u j+1 B jd+l (t) + w j+1 B j+ld+2 (t) - (uj + w j+2 ) B jd+2 (t); (15) 



and 



d_ 
dt 



d 



C jd+1 (t) = itj-iCj-ij+i^) + w j+2 C jd+2 (t) - (u j+1 + Wj) C jd+1 (t) 



u 



J L,i+l(*) 



j t C jtj+2 (t) = u j+1 C jd+l (t) + w j+1 C j+ld+2 (t) - (uj + w j+2 ) C jd+2 (t) 

+ w j+1 B j+lij+2 (t). (16) 

In the limit of t — > 00, again following Derrida's suggestions we introduce the 
ansatz 

Bj,k{t) bj t k, Cj t k(t) — a,j k t — > Tj jk . (17) 

Note that the parameters bj s k, (ij,k and Tj^ are periodic, i.e., bj t k = bj + ^ t k+N, <ij,k = 
dj + N,k+N, and Tj k = Tj +N>k +N- Now define two new functions, 

fj-! = Wjbjj+i ~ Uj-ibj-^j+t, ffl-y = w j+2 b jd+2 - u j+1 b jd+1 . (18) 



At steady state 



dBj, k (t) 
dt 



0, and Eqs. (|15j) transform into 



= u j - 1 b j - 1J+1 + w j+2 b jd+2 - (u j+ i + wj) b jd+1 , 
= Uj +1 b jd+1 + w j+1 b j+ld+2 - (uj + w j+2 ) b jd+2 . 
Substituting (|18j) into these equations, we obtain 



/ 



j - "'./• ~ u j°j,j+2 

(2) 



fo, 



f j+1 - w j+2 b jd+2 - u j+1 b jd+1 - f 



o- 



(19) 



(20) 



where fo is a constant. Then it can be shown that f^ = fj_ x = fo = fj = ff+v This 
leads to the following expression for 

~fo , w j+2 -fo 



r(i) 



bj,j+2 



"I u j,j+2 
Uj+i Uj+i Uj + i 



1 + 



w j+2 



U ; 



-fo . Wj+l- 



-fo 



Uj 



Uj 



J 3+l,j+2 



U; 



Wj+i 

u j+2 



W j+ iW j+2 
UjU j+1 



UjU j+2 



(21) 



Dynamic Properties of Motor Proteins 



8 



Solving these equations by iteration, and using the periodicity and the normalization 
condition, 

JV-1 

we finally derive 

JV-1 

Rn 



j=0 



(22) 



(23) 



where 



U j+ i 



N-lj+k-X 

i+E n 

fe=l i=j 
JV-1 j+fc-1 



fc=i t=j 



+ 



w j+2 



U j+2 



N-lj+k-1 

i+z n 

fc=l i=j 
N-lj+k-1 

i+E n 

fe=i t=j 



To determine the coefficients a.^ and the ansatz (JT7j) is substituted into Eqs. 
(fTKj) in the limit of large times. This yields the following equations 

= Uj- 1 a,j- 1) j +1 + w j+2 a jt j +2 - (u j+1 + w 3 ) a jtj+1 , 

= u j+1 a jJ+1 + w j+1 a j+1J+2 - {uj + w j+2 ) a jtj+2 . (25) 
Also the parameters Tj^ must satisfy 

a j,j+i = + w j+2 T jd+2 - (u j+1 + Wj) T j>j+1 - Uj_ibj_ 1>j+1 , 



- u j+iTj,j+i + tu i+i^i+ij+2 



(Uj + w j+2 ) T jd+2 - w j+1 b j+hj+2 . (26) 

Comparing Eqs. (j25j) with Eqs. (fT9|) we conclude that aj t k = Abj^- The coefficient A 
can be found using the normalization condition IJ22j) and it is equal to 

N-l 1 _ ( Y\ N ~ l HiV 

A = ^ r J,i+2 - WjTjJ+l = N _ V^=0 uj) 

n Rn Rn 

3=0 

To calculate the coefficients T 3 ^ we introduce another set of auxiliary functions 

yf+i = w j+2 T jj+2 - u j+1 T jij+1 , yf_ x = WjTjj+i - Uj-iTj- ljj+ i. 
Then Eqs. (|2fi|) can be rewritten in the following form 

(!) ( 2 ) , l 

a.jj+1 = y j+1 - y)^ + u j - 1 b j - 1J+1 , 

(2) (1) 7 

aj, j+2 = y) ' - - w j+1 b j+1J+2 . 
As shown in [23121], these equations are solved to yield the functions y^ and y^ 



(27) 



(28) 



(29) 



y 



(i) 



.4 



N-l 



-a 7 _i 



,j+l + ( 1 + ^2 d + ■*") + 



i=0 



(2) I, 
S/i = U 3 h 3,j+2 + 



.4 

iV 



'JV-1 

E 

, i=0 



Db 



'j+i+l,j+i+2 + 



a 



2- 



(30) 
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where the coefficients C\ and C 2 are arbitrary constants that cancel in the final 
expression for dispersion. Then the final expressions for Tj^ are found to be 



W j+2 



U; 



'N-\ 

n 

J=o 



W : 



N-l 



.(2) 



U j+1 
j+k-1 



N-l 



yj+i 



j+k-i 



e«&. n 



k=l 



Wj+iWj+2 
UiU i+2 



n 

k=l i=j 



W.j 



'N-l 
J=0 3 



W i+1 W i+3 
U i+1 U i+2 

N-l 



U j+2 



N-l 



j+k-1 



j+k-1 



»? , +E«a n 



k=l i=j 



k=l 
U i+1 U i+3 



i=3 



W i+1 W i+3 
U i+1 U i+2 



(31) 



For simplicity the trailing subunit is chosen as a marker for derivation the explicit 
expressions for the drift velocity and dispersion. It can be shown that the same results 
are obtained if the center of mass is used. The position of this particle at any time is 
given by 

oo N-l 



( X (t)) 



d 
N 



E Eo'+ m ) [ p w+i(«, t) + p j>j+ 2(i, t)\ 



l=—oo j=0 



d 

N 



N-l 



(32) 



3=0 



Utilizing this result along with the Master equations we get 

l=—oo 3=0 

- W 3 P 3,3+l( 1 ' t) + Wj+lPj+ij+2{l, t) - UjP jjj+2 {l, t)] . (33) 

Then the average drift velocity (J2J) has a simple form V = where the function A is 
given by Eq. ()27j1 . The final expression for the velocity is 

{\x=o^Y 

V = d± K (34) 



1 



R 



N 



A similar approach can be used to determine the dispersion, which can be done 
with the help of the following equation, 



d 



2 +oo N-l 



< x2 w> = ^EE^'+ iV/ ) 2 [*wm) + *w>*)i ■ (35) 

l=—oo j=0 

The time evolution of this quantity, again applying the Master equations (fTTj) . is given 
by 

2 {n) E 



d(x 2 {t)) 
dt 



UjC j>j+2 (t) - WjC jtj+1 + ^ (ujB j>j+2 + WjB j>j+1 ] 
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Then after some algebra the expression for dispersion is written as 

2 N-l 



D 



E 

3=0 



~M T j,j+i + Tjj+z)] ■ 



And using the definition © we obtain the final formula: 

d fdU N + VS N N + 2/ 



D 



where 



U, 



N 



N 



N-l 

E 

N-l 



R 2 N 



■V 



A 



(37) 



(3? 



N-l 



j=0 i=0 

JV-1 j'-fe+l 



si(j') 



s 2 (j) 



If ; 



JV-lj'-fc+l 

k=l %=] J 
N-lj-k+1 

WiW i+2 , W j+2 



e n 

fc=i t=j 

JV-l j'-fc+l 



Wi- 2 Wj 



l + V TT -^i±L + ^ ( l + V TT ^i±i I .(39) 

In this mechanism the bounds on the randomness parameter r can be estimated by 
calculating the dynamic properties for the simple limiting case when Uj = u and Wj = 
for all j. Then Eqs. fj23|) . (|52 Jl and (|33 )l yield 



P 2iV T/ W M 

-KAr — , Vinch — d^rrz, 



u 2N 
while for dispersion we get from Eqs. (J58|) and ([3~9*]) 

3N 



27V 2 (iV + l) 

•JAf — ~ , Un 



D 



inch 



U 



(40) 



(41) 



This analysis leads to r = 1/2N which is the smallest possible value for this parameter. 
For any other set of transition rates {uj,Wj} the velocity is always smaller and the 
dispersion is larger, giving the general inequality for the inchworm mechanism 

(42) 



2N 



Although we considered here only the case of m — 2, the method can be extended 
to include the inchworm models where the particles can be found more than 2 sites 
apart. 



2.4- Comparison of Two Mechanisms 

The existence of exact analytical expressions for the dynamic properties of motor 
proteins with two heads in the hand-over-hand and the inchworm mechanisms allows us 
to analyze and compare these mechanisms very efficiently 
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Consider first the simplest N = 1 models. Then the average velocity and dispersion 
for the hand-over-hand mechanism is given by 

V hoh = d(u-w), D hoh = d 2 (u + w)/2, (43) 

The corresponding expressions for the inchworm mechanism can be obtained from Eqs. 

(ESI), 422), (EH), §B and (USD 

V inch = d{u - w)/2, D inch = d 2 {u + w)/8. (44) 

Thus the mean velocity in the inchworm model is only half of the velocity in the hand- 
over-hand mechanism, while the inchworm dispersion is only a quarter of the hand-over- 
hand value. 

More interesting case is iV = 2 models where the average velocity and dispersion 
for the hand-over-hand mechanism are [S2 EDI EI] 

V hoh = d(u Q Ui-w Wi)/a, Dhoh = ^d 2 [(u ui + w wi) - 2(V hoh /d) 2 ] /a, (45) 

where a — Uq + u± + Wq + wi. For the inchworm mechanism the expression for the mean 
velocity can be written as 

T/ . (MpMi) 2 - (w Wi) 2 

Knch = d r r- r, (46) 

2(X(u U 1 + WqWi) + (UqWq - UiW^IUq + W - Ui - Wi) 

while the explicit formula for dispersion is very bulky and it will not be presented here. 
Instead, we analyze the dependence of the dynamics properties of motor proteins on 
external forces using Eqs. ©. 

Force- velocity curves for different mechanisms are presented in Fig. 4. It can be seen 
that the velocity for the inchworm mechanism is always smaller then the corresponding 
curve for the hand-over-hand mechanism, although the stall forces are the same. This 
can be explained by recalling that the stall force is a thermodynamic parameter for the 
sequential chemical kinetic models It is equal to the free energy difference between 
two consecutive binding sites divided over the step size d. Both the free energy difference 
and the step size are the same for the hand-over-hand and the inchworm mechanism, 
and this leads to the same value of the stall force. 

The properties of dispersions for two mechanisms at different external loads, as 
shown in Fig. 5, are similar to the mean velocities. The particles that move via the 
inchworm mechanism fluctuate much less than the motor proteins utilizing the hand- 
over-hand method. This behavior is expected since one of the motor subunits lowers 
the stochastic fluctuations of another motor head in the inchworm mechanism. 

It is also interesting to compare the dimensionless function randomness for each 
mechanism: see Fig. 6. These results suggest that the motor proteins in the inchworm 
mechanism move slower and fluctuates less than the particles in the hand-over-hand 
mechanism, but the relative decrease in the fluctuations is larger than the relative 
lowering of the average speed. This observation, that is correct for any N, is very 
important and it can be used for the experimental discrimination between different 
mechanisms of motor protein's motility. 
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Figure 4. Force- velocity curves for N = 2 model with u = 10 s 1 , iti = 100 s 1 , 
w = 1 s _1 , wi = 10 s - \ 6»+ = 9q = 6»+ = 0f = 0.5 and d = 8.2 nm. The solid 
line describes the hand-over-hand mechanism, while the dashed line corresponds to 
the inchworm mechanism. 



We propose the following procedure to determine the mechanism of the motion 
of motor proteins using only experimental measurements. First, from the independent 
bulk biochemical kinetic experiments determine the number of rate-limiting intermediate 
states. This information provides the size of the period, i.e., the parameter N. Second, 
from the high-precision single-molecule trajectories extract the velocity and dispersion 
for different [ATP] and different external forces. Such data can be obtained from the 
single-molecule optical trap experiments [01 13 IHl M CHI EI]- In the final step, analyze 
the randomness. If for some system this procedure yields r < 1/N, and the known 
number of intermediate states is N, it indicates that the motor protein cannot move 
by the hand-over-hand mechanism. However, for r > 1/N both mechanisms are still 
possible. 

3. Summary and Conclusions 

The dynamics of motor proteins that move along the linear molecular tracks is discussed 
by taking into account the molecular structure and analyzing in detail two possible 
mechanisms of motility. The motor proteins are viewed as two interacting particles that 
correspond to different motor domains in many conventional molecular motors |3|. The 
explicit expressions for the velocity and dispersion are obtained for the hand-over-hand 
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Figure 5. Dispersion as a function of external loads for different mechanisms of the 
motion. The solid line describes the hand- over- hand mechanism, while the dashed line 
corresponds to the inchworm mechanism. The parameters are the same as in Fig. 4. 




Figure 6. Randomness at different external forces. The solid line describes the hand- 
over-hand mechanism, while the dashed line corresponds to the inchworm mechanism. 
The parameters are the same as in Fig. 4. 
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mechanism, when the motor heads pass each other in the alternate fashion, and for the 
inchworm mechanism, when one motor domain is always ahead of another one. 

The exact calculation of the dynamic properties of molecular motors in the hand- 
over-hand mechanism are performed by mapping the dynamics of two particles into the 
one-particle system, for which the dynamic properties are known exactly. It shows that 
in this case the dynamics is identical to the motion of the single free motor domain 
on the same biochemical pathway. The situation is very different for the inchworm 
mechanism. In this case we derived the exact analytic expressions for the velocity and 
dispersion by generalizing the single-particle Derrida's method [2H I2D E3 121] to the 
system with two interacting particles. 

Comparing the dynamics of molecular motors in two different modes, we conclude 
that the proteins in the inchworm mechanism move slower and fluctuate less than the 
particles in the hand-over-hand mechanism. Our results also indicate that the relative 
decrease in dispersion, expressed via the randomness parameter, is smaller for the 
inchworm mechanism. One suggestion is to use this observation for the analysis of 
experimental data on motor proteins. The method of possible discrimination between 
two mechanisms of motor protein's motility based on experimental observations is 
presented. In addition, the effect of external forces on the dynamic properties of 
molecular motors in two mechanisms is also discussed. 

The dynamic properties of motor proteins that move through the inchworm 
mechanism are obtained via the two-particle calculations. However, the average velocity 
and dispersion could also be obtained by mapping the system with two motor domains 
into the system with only one particle, for example, the center of mass of the molecule. In 
general, the inchworm model where the distance between the individual motor domains 
is not larger than m sites, can be mapped into the motion of a single particle onm-1 
parallel biochemical pathways, for which the dynamic properties are known exactly |25j . 

Our analysis of motor protein dynamics is rather very simplified since we considered 
the molecules in which subunits interact only through the hard-core exclusion potential. 
However, the heads in motor proteins coordinate their motion and thus interact much 
stronger then otherwise might be expected It will be interesting to investigate the 
motor proteins with more realistic interactions between the subunits. The theoretical 
method used here seems capable of investigating the more realistic systems of molecular 
motors. 
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